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We investigate a quantum heat engine with a working substance of two particles, one with a spin 
1/2 and the other with an arbitrary spin (spin s), coupled by Heisenberg exchange interaction, and 
subject to an external magnetic field. The engine operates in a quantum Otto cycle. Work harvested 
in the cycle and its efficiency are calculated using quantum thermodynamical definitions. It is found 
that the engine has higher efficiencies at higher spins and can harvest work at higher exchange 
interaction strengths.The role of exchange coupling and spin s on the work output and the thermal 
efficiency is studied in detail. In addition, the engine operation is analyzed from the perspective 
of local work and efficiency. The local work definition is generalized for the global changes and 
the conditions when the global work can be equal or more than the sum of the local works are 
determined. 
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I. INTRODUCTION 


The investigations of heat engines in the quantum 
regime, or quantum thermodynamics, has become an ac¬ 
tive area of research in the last decade fll-[42l| . A quan¬ 
tum heat engine (QHE) uses a quantum working sub¬ 
stance to harvest work in a quantum thermodynamical 
cycle SHI- Three level masers can be considered as 
the first QHEs 1J. Prototype quantum systems, such 
as two level @-@1 and multilevel particles 17-91, coupled 
spins GM3 , and harmonic oscillators 1-1 ll are con¬ 
sidered as quantum working substances. Circuit and cav¬ 
ity quantum electrodynamics systems 351-137}. quantum 
dots [381, quantum Hall edge states 39], cold bosonic 
atoms 40j, optomechanical systems [4l|, and a single 
ion [33:, 42| have been proposed to realize QHEs; while ul¬ 
tracold atoms are porposed for work measurements HE . 
In addition to the studies focusing on the quantum prop¬ 
erties, such as quantum coherence and correlations, of 

the working substance HE, HE HE HE HE HE 1271 - 130} . 

there are explorations of the quantum heat reservoirs as 

weii he mil . 


In the present contribution, we assume classical heat 
reservoirs, and consider two interacting particles, one 
with a spin 1 /2 and the other with an arbitrary spin (spin 
s), as our working medium. The particles are assumed to 
be in an external magnetic field and they interact with 
each other by Heisenberg exchange coupling. The two 
spin 1/2 case of this model has been a s ubj ect of much 
attention in quantum thermodynamics [10T271130} ■ An 
appealing property of the Heisenberg model is that the 
quantum Otto engine efficiency can be enhanced at a crit¬ 
ical exchange interaction between two spins 1/2 HU- We 
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consider the arbitrary spin s as another control parame¬ 
ter next to the exchange coupling and explore its influ¬ 
ence on the performance of the QHE. Such higher spin 
Heisenberg models could be implemented for QHE opera¬ 
tions in nuclear magnetic resonance (NMR) systems 441. 
Among typical quantum thermodynamical cycles [2|, [3| 
we choose to operate our QHE in the Otto cycle as it con¬ 
sists of less demanding processes to implement in com¬ 
parison to other quantum cycles and proposed in various 
systems for implementations HE HE HE HE ■ 

The system we consider can be interpreted as the 
central spin (Gaudin) model with homogenous cou¬ 
plings m The large spin s in our model plays the role 
of a collective spin bath HE consisting of 2s spins 1/2, 
which are homogeneously coupled to a central spin 1/2. 
Recent studies revealed that a spin 1/2 ensemble can be 
used as a heat reservoir to a single spin 1/2, if it consists 
of at least two spins 1/2 HE- In our case there is an ad¬ 
ditional heat reservoir coupled both to the central spin 
and to the collective spin bath. In such a case, the central 
spin 1/2 can always be thermalized while spin s cannot if 
s > 1/2. Even when spin s is not in thermal equilibrium, 
the total system is always fully thermalized. The cou¬ 
pled spin 1/2 and spin s model is hence far from a trivial 
extension of coupled spin 1/2 system but an intriguing 
generalization. 

Consequences of the coupling heat and work reservoirs 
to interacting asymmetric spins lead to surprising results 
which cannot be expected and understood by the knowl¬ 
edge accumulated from the models of coupled spins 1/2. 
For example, we calculated the work output and effi¬ 
ciency of our model QHE. Our results show that with an 
arbitrary spin s, one can extract more work with hi ghe r 
efficiency than the two spins 1/2 case given in Ref 111} . 
Especially the upper bound of efficiency for the two spins 
1/2 given in Ref. 0 can be beaten by an arbitrary spin 
s. In addition, local and global thermodynamics of asym¬ 
metric spins exhibit peculiar differences from those of two 
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spins 1/2. Asymetric two spins with s > 1/2 can act as 
a QHE even in the ultrastrong coupling regime, contrary 
to the two spins 1/2. Furthermore, spin 1/2 can be a re¬ 
frigerant in this regime, while spin s dominates the QHE 
behavior for the total system. Local thermodynamics 
can be explored deeper by the concepts of local heat and 
temperature. Surprisingly local temperature is not appli¬ 
cable to spin s > 1/2 when it is coupled to the spin 1/2. 
The local temperature of spin 1/2 is always well defined. 
It can be controlled by spin s and can be made negative 
or larger than the temperature of the heat reservoir. 

In addition, we introduce generalized definitions of lo¬ 
cal and cooperative work to explore local thermodynam¬ 
ics even when global parameters in a bipartite system 
varies. The earlier definitions are limited to the varia¬ 
tion of local parameters only. Our formalism explicitly 
relates the quantum covariance, or quantum fluctuations 
and quantum coherence, to the work output of QHE 
of any interacting bipartite system. Within the devel¬ 
oped framework, conditions for the violation of the ex¬ 
tensive behavior of the global work as well as a measure of 
quantum cooperativity in terms of covariance naturally 
emerge. While the formalism is model independent, we 
provide an illustration by applying it to the asymmetric 
spin system. 

The paper is organized as follows. In Sec. m we in¬ 
troduce our model QHE. The results for the global and 
local engine operations are given in Secs. Mil and IIV1 re¬ 
spectively. A general discussion on the relation between 
global and local work is given in Sec. El The conclusions 
are stated in Sec. ED 

II. MODEL QUANTUM HEAT ENGINE 

The working substance of our QHE consists of two 
spins in a homogeneous magnetic field, coupled to each 
other with a Heisenberg exchange interaction and it is 
described by a Hamiltonian [lD, 23, 48]: 

H = 8Js a .S b +2B(s z a + S z b ), (1) 

where H = 1 is taken. sa = ( s a; s A; s a)> = 
(S B , S B , S B ), s l A and S B (i = x,y,z) are the spin 1/2 
and spin s operators, respectively. Here, we label the 
spin 1/2 and spin s with A and B , respectively. The 
factor B in the second term of the Hamiltonian denotes 
the external homogeneous magnetic field applied along 
the z axis. We take p,B = 1 and assume there is no or¬ 
bital angular momentum so that the gyromagnetic ratio 
7 is the same for both spins, 7 = 2. J (> 0) is the 
anti-ferromagnetic coupling constant. Here we restrict 
ourselves to s = 1/2,1, 3/2, 2, 5/2, 3. 

The eigenvalues E n of the model Hamiltonian are tab¬ 
ulated in Appendix. In thermal equilibrium with a heat 
bath at temperature T the density matrix p of the work¬ 
ing medium can be written as 

P = Y J Pn\*n){*n\- ( 2 ) 


The occupation probabilities of the eigenstates IT^) 
are P n = exp\—E n /T)/Z (ks = 1) and Z = 
exp (— E n /T ) is the partition function. 

We consider the working medium described by the 
Hamiltonian in Eq. 0 undergoes a quantum Otto cycle 
which consists of two quantum adiabatic and two quan¬ 
tum isochoric processes. The adiabatic branches involve 
the change of magnetic field between two chosen values 
(Hi — > B 2 —S- B\) at a fixed coupling strength, J. The 
details of the cycle are described below. 

Stage 1. This stage is the quantum isochoric process , 
where the working medium with external magnetic field 
B\ and coupling constant J interacts with a heat bath 
at T = T\. The interaction takes long enough, so that 
the working substance falls into a steady state given by 
Eq. (0 with occupation probabilities P n and energy lev¬ 
els E n . Stage 2. The working medium undergoes a quan¬ 
tum adiabatic process , in which the interaction between 
the system and the heat bath is turned off and the mag¬ 
netic field is changed from B 1 to B2. The quantum adia¬ 
batic theorem is considered to hold (provided the process 
is slow enough) [lj| , so that the occupation probabilities 
remain unchanged, while the energy levels change from 
E n to E n due to the change in the magnetic strength. 
Stage 3. This process is almost the reverse of Stage 1, 
where the working medium is in contact with a cold heat 
bath at T = T 2 (Tf > X 2 ). Reaching equilibrium with the 
bath changes the energy probabilities to P n with B = B 2 , 
T = T 2 and J in Eq. 0. Stage /. The system undergoes 
another quantum adiabatic process with changing B2 to 
B\ ( E n to E n ), while keeping P n the same. 

From the generalization of the first law of thermody¬ 
namics to quantum mechanical systems @-S , the heat 
exchanges in Stages 1 and 3 are, respectively, given as 

Ql = y] E n ^ P n ~ Pr^j , 

n 

Q 2 = J2 E 'n( P n-Pn)- (3) 

n 

The work is performed only in the adiabatic branches 
of the quantum Otto cycle. Due to the conservation of 
energy, the net work done by the QHE can be written as: 

W = Qi + Q 2 

= E - K) ( P n - Pn) , (4) 

n 

where W > 0 signifies the work performed by the QHE 
with operational efficiency r) = W/Q\. To harvest posi¬ 
tive work by the engine, we consider Qi > —Q 2 > 0 to 
conform to the second law of thermodynamics. 

By using the tabulated eigenvalues E n of H in Ap¬ 
pendix and the probabilities given by the thermal oc¬ 
cupation numbers in Eq. ©, the work output and the 
efficiency of the engine can be calculated analytically. 
The analytical expressions are not very illuminating and 
will not be displayed here for brevity. We call the work 
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done by the engine given by Eq. © and its efficiency rj 
as the global work and global efficiency, respectively, to 
distinguish them from the local work and efficiency of 
individual spins, described later in the text. 


III. GLOBAL WORK AND EFFICIENCY 


Before presenting our results, we would like to review 
some of the main results in Ref. (ll| where the authors 
investigated the same Hamiltonian in Eq. m but for two 
spins 1/2. The conditions in which the coupled engine ef¬ 
ficiency can be higher than the uncoupled one have been 
determined. Specifically, an upper bound rjb to the effi¬ 
ciency 77 of the quantum Otto engine has been obtained 
as 


7? < Tjb 


1 ~ B 2 /Bi 
1 - 4 J/B 1 


< Vc 


(5) 


where the upper bound is always less than the classical 
Carnot efficiency (r; c = 1 — T2/T1). 

In Fig. [□ we investigate the role of spin s on the 
performance of the coupled quantum Otto engine. We 
plot the global work in Fig. [Qa) and global efficiency in 
Fig- Hlb), as a function of exchange coupling strength J 
for Bi > i ?2 and s = 1/2, 1, 3/2, 2, 5/2,3. For the uncou¬ 
pled engine (J = 0), the engine efficiency can be calcu¬ 
lated as r)j—o = 1 — B2IB\ which is independent of spin 
s as can be seen in Fig. QJb). The coupled engine per¬ 
formance can be higher than the uncoupled one; both W 
and rj first increase to certain maximums as a function of 
J and then drop to zero. The role of spin s on the global 
work and efficiency is found to shift the maximums and 
the positive work conditions (PWCs) to the weak cou¬ 
pling regimes; accordingly the coupled Otto engine with 
high spin s can produce higher work with higher efficiency 
than the lower spin s, below a certain sufficiently weak 
coupling strength (for instance, J <« 0.12 in Fig. [TJ. 
Especially, the engine with s > 1/2 can violate 775 as 
indicated by the dashed line in Fig. [Tf b). 

The mutual relationship between the work output and 
efficiency is demonstrated by the characteristic curve in 
Fig. [2l for the same magnetic field and temperature val¬ 
ues, and for the same coupling strength range as in Fig.Q] 
It can be deduced from Fig. [2] that the efficiency at max¬ 
imum work output as well as the work at maximum effi¬ 
ciency are not notably affected by the spin s of the work¬ 
ing substance. It seems that the higher spin s leads to 
higher efficiency and work output at the weak coupling 
regime. We should stress here that this is not the general 
conclusion; for differently tailored parameters, the max¬ 
imum of work output and the efficiency can slightly be 
influenced by the spin s. 

In Fig. [1} we have restricted ourselves to the weak cou¬ 
pling regime, specifically J £ [0,0.5], and now we focus 
on the strong coupling region. It is possible to show that 
beyond this limit, i.e., J > 0.5, the working substance of 
two spins 1/2 cannot do positive work, since it violates 




FIG. 1: (Color online.) Dependence of global work W (a) 
and efficiency rj (b) on coupling strength J for temperatures 
Ti — 1, 7■> = 0.5, and magnetic fields B\ = 4, B 2 = 3, 
and spins s = 1/2 (black line), and s = 1 (red line), s = 
3/2 (blue line), s = 2 (green line), s = 5/2 (yellow line) and 
s = 3 (magenta line). The dashed line in (b) indicates the 
upper bound r/i, of the global efficiency given in Eq. ([5]) for 
the case of spin 1/2 pair. For the above parameters, we have 
r) c = 1 - T 2 /T 1 = 0.5 and r, J=0 = 1 - B 2 /Si = 0.25. All 
quantities plotted are dimensionless. In all figures, we use a 
unit system where h = 1,pb = 1, fcs = 1 and use Ti as our 
scaling parameter. 


the PWC given in Ref. EH- It is reasonable to assume 
that the change of energy gaps in the adiabatic stages by 
the change of magnetic field cannot contribute in the di¬ 
rection of total positive work gradient when J > 0.5. On 
the other hand, for the case of pairing spin 1/2 and spin s 
with s > 1 / 2 , the role of energy gaps in the work extrac¬ 
tion can be dramatically changed after a critical value 
of J and the engine can reproduce useful work. This is 
shown in Fig. [3] where the global work and efficiency are 
plotted as a function of J up to the very strong cou¬ 
plings. As shown in Fig. [3j the positive work re-emerges 
after a critical value of coupling strength. Increasing the 
spin s value shifts the critical J towards the weak cou¬ 
pling regime. The efficiency and the work output are less 
in the strong coupling regime. Since the corresponding 
thermodynamical quantities are invariant under uniform 
energy shifts [§], the coupled spin 1/2 and spin s model 
in the limit of very large coupling strengths (i.e., J — 00 ) 
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FIG. 2: (Color online.) Mutual relation of global work W 
and efficiency r/ for the same magnetic field and temperature 
parameters and the coupling range as in Fig. |T] The curves 
for each spin s nearly coincide. 




J 


FIG. 3: (Color online.) The work output (a) and efficiency 
(b) in the broader range of J, including the strong coupling 
region for the same parameters and spin s as in Fig. [T] The 
direction of arrow in (b) indicates the lines in the order of 
increasing spin s from s = 1 to s = 3. The curves are the 
same with those in Fig. [l]in the weak coupling regime. Note 
that after J « 0.5, the engine cannot produce positive work 
for the case of coupled spins 1 / 2 . 


can be mapped into a multilevel system with energy lev¬ 


els { 0 , 2 B, 4 B, ..., (2s — 1)2 B} where 77 = 0 for s = 1/2, 
while 77 = 1 — B 2 /B 1 for s > 1/2. This explains the be¬ 
havior of the efficiency in Fig. [31(b) where 77 converges to 
the spin independent value of 77 = 1 — B 2 /Bi for s > 1/2 
and 77 = 0 for s = 1/2 in the deep strong coupling regime. 


IV. LOCAL WORK AND EFFICIENCY 

In this section, we investigate how the spin 1/2 and 
spin s individually undergo the engine operation. This 
can be done by the analysis of local heat exchanges be¬ 
tween the local spin and the reservoir (Tlj • The local heat 
exchanges in the isochoric branches of the Otto cycle can 
be expressed as the change in the local density matrix for 
a given local Hamiltonian. Let q\ (q 2 ) : with i = A,B, be 
the local heat transferred between the ith spin and the 
hot (cold) heat bath. Then the explicit expression of q\ 
(q 2 ) reads as 0 : 

q\ = Tr[{ Pi ~ Pi )Hi], 

q l 2 = Tr[{ Pi — (6) 

where pt ( Pi ) is the reduced density matrix for the ith 
spin at the end of stage 1 (3) and Hi ( H t ) is the local 
Hamiltonian during the first (second) isochoric process. 
The local Hamiltonians can be written as Ha = 2 Bs z A 
and Hg = 2 BS b for the spin 1/2 and spin s, respectively. 
The local work done by the ith spin is then written as 
Wi = q\ + q 2 . 

The local works vja and wb, done by the spin 1/2 and 
spin s, respectively, are plotted as a function of coupling 
strength J in Fig. [4] for s = 1/2,1, 3/2,2, 5/2, 3. The 
analytical calculation of the global and local works yields 
that W = wa + wb, the total work is the sum of local 
efforts. For further insight, it is possible to calculate 
the relation between the global and local heat exchanges, 
which is found to be 

Qi = q\ + qf +8JV S , 

Q 2 = q A + gf - 8 JV S1 (7) 

where V s = Tr[(p — p )sa-Sb ], p (p ) being the global 
thermal density matrix at the end of stage 1 (3) given 
by Eq. ©• The term V Sl given by expectation value of 
interacting part, is related to the probabilities of certain 
energy levels at the end of stages 1 and 3. Its explicit 
expression depends on the spin-s but not written here 
explicitly for brevity. The relations in Eq. 0 suggest 
that only the local heat exchange is converted into total 
work output of the Otto cycle, as the last terms in Q\ 
and Q 2 expressions reflect the collective heat intake and 
release which cancel each other. This is consistent with 
the extensive property of the work output of the cycle. 
We should stress here that same conclusion is reached 
for the case of spin 1/2 0 and spin 3/2 pairs 0. The 
extensive property is not a fundamental character of the 
work output and is not always true. Similar analysis in 
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FIG. 4: (Color online.) The local work done by the spin 1/2 
(a) and spin s (b) versus coupling strength J for values T\ = 1, 
T 2 = 0.5, B\ = 4, B 2 = 3, and s = 1/2 (black line), s = 1 (red 
line), s = 3/2 (blue line), s = 2 (green line), s = 5/2 (yellow 
line) and s = 3 (magenta line). 


different conditions reveal that sum of the local works 
is not always equal to the global work [kJ [2f| . We will 
present a more general discussion in the following section. 

For two coupled spin 1/2 case, we have wa = wb since 
PA = Pb and H A = H B a Moreover, for J = 0, wa 
is independent of spin s value, wb depends on spin s for 
J = 0, but this dependence is weak to be visible in the 
scale of Fig. 01 On the other hand, these results are dra¬ 
matically changed when s > 1/2 and J / 0. As shown 
in Fig. 0Ja), wa depends strongly on the spin s. In the 
region J < 0.5, increasing s shifts the PWCs and max- 
imums of wa and wb to the weak coupling regions and 
increases (decreases) the maximums of wa (wb). The 
comparison of local works of both spins shows that, ex¬ 
cept a negligibly tiny range of J, we have wa > wb, that 
is spin 1/2 does more work than the spin s. On the other 
hand, if we change our attention to the strong coupling 
regime where J > 0.5 (Fig.0, this situation is completely 
reversed; as shown in Fig. 0a), wa < 0 for each spin s, 
while wb can be non-zero for s > 1/2 (Fig. [5](b) ). From 
an analytical calculation of global and local works in the 
deep strong coupling regime (i.e., J —>• 00), it is possible 
to show that W = —(2s + 1 )wa = (2s + l)/(2s + 2 )wb- 
This indicates that spin s is solely responsible for the 
realization of our QHE in the strong coupling regime, 


where wb > 0 and wa < 0 in the regions W > 0. 




FIG. 5: (Color online.) The local work done by the spin 1/2 
(a) and spin s (b) in the strong coupling region where J > 0.5 
for the same parameters and spin s as in Fig. [4] The curves 
are the same with those in Fig.0]in the weak coupling regime. 
The direction of arrow in (b) indicates the lines in the order 
of increasing spin s from s = 1 to s = 3. Note that after 
J fts 0.5, wa = wb < 0 for s = 1/2 as given black line in (a). 

For further insight of the local performance, we ana¬ 
lyze the local efficiency and the local heat flow of spin 
1/2 and spin s. From the local work and heat ex¬ 
change definitions in Eq. ©, for the considered case (i.e., 
Bi > B 2 > 0 IH), we always have q\ > — q l 2 > 0 when 
Wi > 0; that is local heat always flow in the direction of 
global heat gradient when the local work is positive. For 
two spins 1/2, local spins are always heat engines when 
the global work W > 0, since wa = wb = w, W = 2w 
and W 1 w > 0. On the other hand, this is not necessarily 
true for the two coupled asymmetric spins. One of the 
spins can be refrigerant when total system operates as 
a QHE. From the comparison of Fig. [T] and Fig. 01 one 
can deduce that there is a critical J depending on s up 
to which the local spins are heat engines. Beyond the 
critical coupling strength, as can be seen in Fig. 01 we 
have wb > 0 and wa < 0 (g) 4 < 0 and qQ > 0) when 
W > 0. Here spin 1/2 acts in the opposite direction of 
global work gradient and it is a refrigerator, although 
total system and the spin s > 1/2 are the heat engines. 

For the local heat exchanges in Eq. ©, we have the re¬ 
lation q\ = —{Bi/B 2 )q 2 : so we found that the individual 
spins undergo the heat cycle with the same and constant 
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local efficiency: pa = Vb = Wi/q\ = 1 — B 2 /Bi, which is 
independent of spin s and equals to the global uncoupled 
engine efficiency. In the regime where spin 1/2 operates 
locally in the refrigeration cycle it is more appropriate to 
consider local coefficient of performance, which is also a 
constant €a = q 2 /\wa\ = B 2 j(B\ — B 2 ). 




FIG. 6: (Color online.) The local temperature of spin 1/2 
versus the parameter s at J = 0.1 (a) and at J = 4.0 (b). Here 
Ta ( Ta ) denotes the local temperature for the hot (cold) heat 
bath cases with Ti = 1.0 (T 2 = 0.5) and Bi = 4 (B 2 = 3). 


Now we focus on the temperatures of local spins at the 
end of two thermalization stages and investigate the role 
of coupling J and spin s on the local temperatures. By 
taking two different energy levels, Ei, and their proba¬ 
bilities, Pi (obtained from the reduced density matrix in 
Eq. ©), one can define an effective temperature for the 
local spins A and B as: 


T k 


Ei — Ej 

In Pj — In Pi ’ 


k = A,B. 


( 8 ) 


In the absence of coupling between spins (i.e., J = 0), 
the local temperature of both spins are equal to the heat 
bath temperature, since we have p = pa ® Pb, where 
Pi (i = A,B) are the local thermal density matrices at 
the heat bath temperature. On the other hand, this sep¬ 
aration would not be possible in the case of nonzero J 
which can make the spins to locally thermalize to differ¬ 
ent temperatures than the heat bath (llj| . Analysis of 
the local temperature of spin s > 1/2 reveals that the 


reduced state of spin s is a non-equilibrium steady state; 
where there is no unique effective temperature applica¬ 
ble to all the energy level pairs [2j. For the two-level 
system (spin 1/2) it would always be possible to define 
an effective temperature. 

Interplay of interactions and multilevel nature of large 
spin on the thermalization can be explained at more in¬ 
tuitive level as follows. Our model can be interpreted 
as a central spin coupled to an ensemble of 2s spins 
1/2 [45], [46}. Both the collective spin bath and the cen¬ 
tral spin are also coupled to a heat reservoir. The heat 
bath and the spin s > 1/2 can always fully thermalize 
the spin 1/2; while spin s > 1/2 cannot be thermalized 
when it is coupled to the spin 1/2. This intuitive rea¬ 
soning based upon the size difference of subsystems is in 
parallel with the conclusion in Ref. (47) which states that 
the minimum number of spins 1 /2 required to thermalize 
a single spin 1/2 is two. 

The effective temperature Ta of spin 1/2 with respect 
to spin s at the end of first ( Ta) and second (Ta) 
isochoric processes for weak and ultrastrong coupling 
regimes are shown in Fig. (HKa) and in Fig. Ob), respec¬ 
tively. At weak coupling regime (J = 0.1), where the spin 
1/2 is a heat engine as shown in Fig. [4ja) , Ta is always 
larger than the heat bath temperature. Spin s, acting as 
an additional reservoir next to the heat bath, can heat the 
central spin 1/2 to higher temperatures. For two spins 
1/2, increasing J makes the occupation probabilities of 
energy levels of reduced spins equal so that effective tem¬ 
perature becomes infinitely high at the ultrastrong cou¬ 
pling regime. On the other hand, for s > 1, after a criti¬ 
cal coupling strength the population in high energy level 
(£q = B) exceeds the population in lower energy level 
(Ei = —B ) so that local temperature of the spin 1/2 
becomes negative as shown in Fig. Ob). After the first 
and the second adiabatic stages, the local temperatures 
of the spin 1/2 change to ( B 2 /B{)Ta and ( Bi/B 2 )Ta , 
respectively. Fig. O a ) indicates that T% > (Bi/B 2 )Tj( 
and T^ < (B 2 /B\)Ta , so that g) 4 > 0 and q 2 < 0; hence 
wa > 0 as discussed above. In the deep strong coupling 
regime, Fig. Ob) elucidates that T% < (Bi/B 2 )Ta and 
T^ > (B 2 /Bi)Ta , so that gj 4 < 0 , q 2 > 0; consequently 
wa < 0, so that spin 1/2 is a refrigerator. 


V. GENERAL RELATIONS BETWEEN 
GLOBAL AND LOCAL WORK 


We have seen in Sec. Ev] that the global work has 
an extensive property and can be written as a sum of 
the local works done by the individual spins. This con¬ 
clusion strictly depends on the paths, or the methods, 
we choose to operate the engine cycle. In the adiabatic 
stages of the quantum Otto cycle, we varied the homoge¬ 
neous magnetic field acting on the spins. We can make a 
general statement that it is not possible break the exten¬ 
sive property of work output of a QHE by only making 
local changes in the adiabatic stages of the engine cy- 
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cle. This simple fact can be quickly proven for a general 
Hamiltonian of a system of a collection of local subsys¬ 
tems, described in the form H = J2H\ oc + -Hint, where 
the non-interacting (local) and interacting (global) terms 
are denoted by H\ oc and H lnt , respectively. The internal 
energy, U = (H) = Tr (pH), of the system with density 
matrix p changes as dU = Tr(pdH)+Ti{Hdp), where the 
first term can be defined as the work done on the system 
and denoted by dW := Tr (pdH). In a strictly quan¬ 
tum adiabatic process we have dp = 0. Accordingly, if 
dH lnt = 0, the global work becomes extensive in terms of 
local works done by subsystems such that dW = ^dw i oc , 
withdu/ioc := Tii oc (p\ oc dHi oc ), where pi oc is the reduced 
density matrix of a particular subsystem found by trac¬ 
ing out the degrees of freedom of the other subsystems 
from the density matrix p of the whole system. While 
the global work is extensive under local changes, it can 
still be optimized by the interactions between the subsys¬ 
tems, through the interaction dependence of the reduced 
density matrices p\ oc , which is illustrated by our analysis 
in Sec. Ell and Sec. HU 

Let us now consider a more general situation where 
both the magnetic field and the exchange interaction be¬ 
tween the spins could change. In such a case, Eq. m 
directly shows that the extensive behavior of the global 
work is violated by the simultaneous change of magnetic 
field strength {B\ —> B 2 —> B\) and the exchange cou¬ 
pling strength (Ji —>• J 2 —» J\) in the adiabatic stages 
such that 

W = w A + w B + 8 ( J\ - Ji)T s , (9) 

where V s is defined in Eq. 0 . 

A curious result of Eq. 0 is that when B\ = B 2 
and J\ ^ J 2 , the system can harvest positive work in 
a purely collective manner, as no local work can be done 
by the local systems in constant magnetic field. Since 
there is no change in local Hamiltonians, the total lo¬ 
cal heat exchange is zero. If we take the ratio W/w\ oc , 
where w\ oc = wa + w b is the total local work, as a figure 
of merit measuring the cooperativity in work extraction, 
it is infinite. On the other hand, we can still consider 
a possible generalization of the local work definition in 
Ref. [Tlj to scrutiny them in a purely interacting cycle 
without explicit local variations. We suggest that a mean 
field Hamiltonian can always be introduced to describe a 
local Hamiltonian of a subsystem. 

To make our discussion concrete let us take a pair¬ 
wise interaction Hamiltonian of the form H = gAB , 
where A and B are operators for two subsystems, and 
g is their coupling constant. The work done on the sys¬ 
tem in an adiabatic stage by the dg variation of the cou¬ 
pling constant can be written as dW = dg(AB ) 7 where 
(AB) = Tr (pAB). If we use mean field Hamiltonians 
Ha = g(B)A/ 2 and Hb = g{A)B/2 for the local Hamil¬ 
tonians then the corresponding local work contributions 
become wa = w B = g(A)(B)/2. Accordingly the global 
work can be expressed as dW = dwA + dw b + dw coop , 
where we introduced a cooperative work term<?u; C oop := 


dg(A,B). Here, the notation (A, B) := {AB) — (A)(B) 
stands for the covariance of A and B as a measure of 
correlations between the subsystems. The net work done 
in the cycle then becomes 

W = wa + w B + w coop , (10) 

where the local and cooperative works are given by 

w A = w B = ^(si - g 2 ){{A) 1 (B } 1 - {A) 2 {B) 2 ), 

Wcoop = (ffi - 92 ){{A, B)i - (A,B) 2 ). (11) 

Here g\ and g 2 are the coupling constants at the end 
points of the adiabatic stages, and the expectation values 
(X)i = Tr (pxX) are evaluated with the reduced density 
matrix p l x of the subsystem X = A, B in the adiabatic 
stage labeled by i = 1,2. With this generalized definition 
of the local work, the cooperativity of the work extraction 
can be characterized by the ratio 

W (A,B ) 1 - {A,B ) 2 

Wioc + {A) 1 {B) 1 -{A) 2 {B ) 2 1 J 

Applying the generalized local work formalism to our 
Heisenberg exchange model QHE, we find the local 
Hamiltonians 

Ha = 2Bsa + 2 8 J s a-(Sb), 

H b = 2BS z b + ^8J{sa).S b , (13) 

which gives the relation between global and local works 
as dW =dwA +dwB +dw c00 p, where 

dw A = 2dB{s z A )+ ^8dJ{s A )-{S B ), 

dw B = 2dB{S z B ) + ^ 8 dJ{s A ).{S B ), (14) 

and(fu> coop = 8 dJ(s A , S B )- From this result we conclude 
that the extensive property of the global work can be 
violated by changing the interaction parameter in the 
adiabatic stages, if the covariance of the interacting spins 
changes as well. If the covariance remains the same, then 
the global work can be expressed as the sum of effective 
local works of the individual spins under the mean field 
description. 

VI. CONCLUSIONS 

We consider a pair of spin 1/2 and spin s coupled via 
Heisenberg exchange interaction under a homogeneous 
magnetic field as the working medium of a quantum Otto 
engine. The influence of exchange coupling and spin s on 
the work output and efficiency of the quantum Otto en¬ 
gine is investigated in detail. The global engine operation 
is also analyzed in comparison to local work contributions 




of the individual spins. It is found that increasing spin 
s at a certain exchange coupling strength can make the 
QHE to produce more work with higher efficiency, which 
can violate the upper bound of efficiency for two coupled 
spin 1/2 particles & Moreover, spin s makes it pos¬ 
sible to realize the QHE at the strong coupling regimes. 
Furthermore, we show that due to the coupling of asym¬ 
metric spins, one of the spins can operate as a refrigerator 
even when global cycle is a heat engine. From the local 
work analysis, it is found that global work is equal to the 
sum of the local works by the individual spins. Although 
in local realm, the spin 1/2 and spin s operate with the 
same efficiency, their local works are found to be signif¬ 
icantly influenced by the spin s. The local temperature 
of spin 1/2 is found to be controlled by spin s and can 
be negative or larger than the temperature of the heat 
baths in the case of non-zero coupling. Finally we dis¬ 
cussed the conditions for the violation of the extensive 
behavior of the global work. We developed a formalism, 
applicable to any coupled bipartite system, generalizing 
the local and cooperative work definitions to the case 
where global changes can be performed in the engine cy¬ 
cle. The general conditions for which the global work is 
not equal to the sum of the local works are given in terms 
of the covariance of the subsystems. 
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Appendix: The Eigenvalues of the Working Medium 

Here we report the eigenvalues of the Hamiltonian (JT]) 
for s = 1/2,1,3/2, 2, 5/2, 3. The corresponding orthonor¬ 
mal eigenstates can also be calculated. We should stress 
here that the eigenstates are system parameter (i.e., J 
and B) independent. Since the discussion of text does 
not require the explicit form of the eigenstates, we do 
not report them here for brevity. 

The eigenvalues for (i, s) system with s = 1/2 are & 
{-6 J, 2 J - 2B, 2 J, 2 J + 2 B}. 

The eigenvalues for (i,s) system with s = 1 are: 
{—B — 8J,B- 8 J, —3 B + 4J, -B + 4 J,B + 4J, 3 B + 4J}. 

The eigenvalues for (|,s) system with s = 3/2 are: 
{—2H-10 J, -10 J, 2B—10J, -2B+6J, -4H+6 J, 6 J, 2B+ 
6J,45 + 6J}. 

The eigenvalues for (i,s) system with s = 2 are: 
{-3B-12 J, -5-12 J, 5—12 J, 35-12 J, -55+8 J, -3 B+ 
8 J, —5 + 8 J, 5 + 8 J, 3 B + 8 J, 5 B + 8 J}. 

The eigenvalues for (|,s) system with s = 5/2 are: 
{-45-14J, -25-14J, -14J, 25-14J,45-14J, -2 B + 
10 J, -45 + 10J, -65 + 10J, 10 J, 2B +10 J, 45 +10 J, 65 + 
10 J}. 

The eigenvalues for (i,s) system with s = 3 are: 
{-55-16 J, -35-16 J, -5-16 J, 5-16 J, 35-16 J, 55- 
16 J, -3 B + 12 J, 3 B + 12 J, -7 B + 12 J, -5 B + 12 J, —B + 
12 J, B + 12 J, 5 B + 12 J, 7 B + 12 J}. 
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